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ABSTRACT 


Results  from  experimental  studies  concerning  the 
natural  frequencies  and  mode  shapes  of  the  thin  truncated 
conical  shell  with  free  edge  conditions  have  not  been  con¬ 
sistent.  Mode  shapes  from  one  study  were  found  to  have  an 
equal  number  of  nodes  at  the  two  free  cone  edges  when 
vibrating  in  a  given  mode  while  in  contrast  a  second  study 
revealed  an  unequal  number  of  nodes  at  the  edges.  The 
theoretical  results  presented  in  this  investigation  are 
obtained  through  the  use  of  a  modified  Galerkin  procedure 
and  are  in  agreement  with  the  equal  node  experimental  data. 
At  the  same  time,  the  theoretical  results  have  been  used 
as  a  basis  for  explaining  the  opposing  experimental  data. 


vii 
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CHAPTER  I 

INTRODUCTION 

The  subject  of  conical  shell  vibrations  has  received 
very  limited  attention  in  the  literature  until  recently. 
Although  the  conical  shell  has  probably  the  simplest  geome¬ 
try,  with  the  exception  of  cylindrical  and  spherical  shells, 
the  mathematical  analysis  of  a  conical  shell  is  more  diffi¬ 
cult  due  to  the  fact  that  the  governing  diffeiential  equa¬ 
tions  of  motion  have  variable  coefficients.  Consequently, 
most  analyses  of  the  vibrations  of  conical  shells  are 
accomplished  through  the  use  of  approximate  methods  and 
numerical  procedures  for  which  large  digital  computers  are 
a  tremendous  aid. 

A  comprehensive  summary  of  significant  previous 
works  is  given  by  Hu  (4)1.  The  majority  of  analyses  con¬ 
cerning  the  natural  frequencies  and  mode  shapes  of  trun¬ 
cated  conical  shells  concern  those  cases  where  the  shell 
is  fixed  at  one  edge  and  free  at  the  other,  or  simply  sup¬ 
ported  at  both  edges.  A  recent  work  by  Platus  (9)  is  an 
example  of  a  modern  treatment  of  the  fixed-free  case, 

1.  Numbers  in  parentheses  refer  to  the  list  of 
references . 
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wherein  a  digital  computer  and  matrix  methods  are  utilized. 
However,  the  case  where  both  edges  of  the  truncated  conical 
shell  are  free  of  external  support  has  not  received  similar 
analytical  attention.  Perhaps  the  reason  is  that  it  is 
difficult  to  choose  functions  which  satisfy  force  type 
boundary  conditions  of  the  free-free  case. 

Recent  experimental  studies  for  the  case  of  free- 
free  edge  conditions  have  been  made  by  Watkins  and  Clary 
(12)  and  Hu,  Gormley  and  Lindholm  (6).  Watkins  and  Clary 
observed  that  at  certain  natural  frequencies  a  greater 
number  of  circumferential  nodes  occurred  at  the  major  end 
than  at  the  minor  end  of  the  shell.  This  observation  has 
generated  some  discussion  in  the  literature  [see  Hu  (5)] 
concerning  the  validity  of  the  results.  In  the  experi¬ 
mental  work  of  Hu,  at  al. ,  the  phenomenon  of  different 
numbers  of  circumferential  waves  at  the  cone  edges  was  not 
observed. 

At  the  time  this  study  was  begun,  the  results  of 
Hu,  at  al.  were  not  available.  The  objective  of  the 
Investigation  at  that  time  was  to  investigate  the  possi¬ 
bility  that  pure  modes  existed  having  nodal  lines  which 
are  not  generators  of  the  cone,  i.e.,  correlate  theoretical 
analysis  with  the  experimental  data  of  Watkins  and  Clary. 
However,  the  analysis,  wherein  the  assumed  displacement 
shapes  permitted  an  unequal  number  of  circumferential  waves 
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at  the  two  edges  of  the  cone,  revealed  an  uncoupling  such 
that  a  pure  mode  having  an  unequal  number  of  nodes  at  the 
cone  edges  is  not  possible. 

With  the  advent  of  the  experimental  results  of  Hu, 
Gormley  and  Lindholm,  a  controversy  arose  since  their  work 
did  not  reveal  different  numbers  of  circumferential  waves 
at  the  two  cone  edges  in  a  given  vibration  mode.  The 
objective  of  the  study  at  hand  was  therefore  altered  toward 
the  direction  of  attempting  to  explain  the  apparent  con¬ 
troversy  between  the  two  experimental  studies. 

In  this  investigation  an  approximate  solution  for 
the  free-free  case  is  obtained  utilizing  a  modified  Galer- 
kin  procedure  wherein  functions  assumed  in  the  solution 
need  satisfy  only  the  displacement  boundary  conditions,  a 
feature  of  considerable  advantage  for  the  free-free  case. 
The  method  accounts  for  errors  in  the  forces  and  moments 
which  result  when  the  assumed  functions  do  not  satisfy 
the  force  type  boundary  conditions. 

The  theoretical  analyses  proceeded  with  the  com¬ 
puter  programming  of  the  computational  procedure  described 
in  Chapter  IV.  Results  from  the  program  for  the  experi¬ 
mental  case  of  Hu,  et_  al.  showed  good  agreement  with  the 
experimental  data  both  in  frequencies  and  mode  shapes. 

This  correlation  of  results  is  presented  in  Chapter  V. 
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A  plausible  explanation  of  the  experimental  results 
of  Watkins  and  Clary,  basea  upon  specific  computer  results, 
is  presented  in  Chapter  VI. 


CHAPTER  II 


EQUATIONS  OF  MOTION  AND  BOUNDARY  CONDITIONS 

-The  differential  equations  of  motion  utilized  in 
this  Investigation  are  those  developed  by  Seide  (10)  from 
energy  principles.  They  are  identical  to  those  used  by 
Uhao-tsien  (2).  Equivalent  equations,  based  upon  elemental 
equilibrium  considerations,  can  be  obtained  using  Flugge's 
(3)  simplified  version  of  the  elastic  lai*  as  applied  to  a 
thin  conical  shell  of  constant  thickness.  The  equations 
apply  to  homogeneous  isotropic  conical  shells  of  constant 
wall  thickness  and  are  based  upon  linear,  small  displace¬ 
ment  shell  theory  and  the  assumptions  of  a  large  mean  diam¬ 
eter  to  shell  thickness  ratio  and  negligible  deformations 
due  to  transverse  shear.  For  "thin"  shell  theory  it  is 
assumed  here  that  h /  (s  sin  a)  «  1,  where  h,  s  and  a  are 
as  defined  in  Figure  2.1.  More  complete  Flugge  equations 
are  presented  in  Appendix  A.  Basically,  the  additional 
terms  found  in  the  more  general  equations  represent  higher 
order  for  the  very  thin  truncated  shells  analyzed  in 

this  study. 

Development  of  the  equations  based  upon  energy 
principles  will  be  presented  here.  Expressions  for  the 


(2.7) 


1  dw  .  1  ^ 

X8=^+t^t;^ 


V  -  1  S  (1  ^w\ 

*s0  sin  a  cFs  's  cJ0' 
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(2.8) 

(2-9) 


Figure  2.2  shows  the  force  and  moment  stress 
resultants  acting  on  a  differential  shell  element.  Stress 
resultants,  Ng,  Ng0,  and  Qg  are  forces  per  unit  of  circum¬ 
ferential  length;  N0,  N0g,  and  Q,0  are  forces  per  unit  of 
longitudinal  length;  Mg  ard  Mg0  are  moments  per  unit  of 
circumferential  length,  and  MQ  and  M0g  are  moments  per  unit 
of  longitudinal  length.  The  equations  relating  the  stress 
resultants  to  strains  and  curvatures  are  as  given  by  Seide 
(10): 

N  =  U  •+  vO  (2.10) 


NG  “ 


-  Dh 


-  Dh 


■'"s  •*  V€9> 

(2.10) 

>e  +  VEs > 

(2.11) 

Eh  .. 

2 ( 1-fv )  Ys0 

(2.12) 

(Xs  +  vx0) 

(2-13) 

(X9  +  vXs) 

(2.14) 

-  -  Mes  =  W 


(2.15) 
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where  E  is  Young's  modulus,  v  is  Poisson's  ratio,  h  is  the 
shell  thickness,  and  D  is  defined  as 

d  =  __Hi  (2.16) 

(l-vd) 

Equations  of  equilibrium  and  boundary  conditions 
can  be  derived  from  energy  principles.  The  total  energy 
of  a  conical  shell,  i.e.,  the  strain  energy  minus  the  work 
due  to  external  forces,  can  be  expressed  for  the  complete 
circular  shell  as 

2tt  s2 

u  =  I  I  '  t|  ( vs  +  M0ce  +  iy  Yse  -  Ms  xs 

r>  c 


where  Fu>  Fv  and  Fw  are  external  forces  per  unit  of  surface 
area  and  the  barred  stress  resultants  are  forces  or  moments 
per  unit  length  on  the  cone  boundaries.  Terms  to  be  inte¬ 
grated  along  a  longitudinal  edge  from  s^  to  Sg  do  not 


11 

appear  since  the  development  here  is  limited  to  a  complete 
shell  8  =  0  to  2t r.  The  terms,  Fu,  Fy  and  F  ,  will  later  be 
replaced  by  inertia  forces  using  D'Alembert's  principle. 

The  energy  expression  in  terms  of  the  displace¬ 
ments,  u,  v,  and  w,  is  obtained  by  substituting  the  equa¬ 
tions  for  stress  resultants  in  terms  of  strains  and  curva¬ 
tures,  Eqs.  (2.10)  to  (2.15),  and  strains  and  curvatures  in 
terms  of  displacements,  Eqs.  (2.1)  to  (2.3)  and  (2.7)  to 
(2,9),  into  the  energy  expression,  Eq.  (2.17).  When  the 
first  variation  of  the  resulting  expression  is  taken  and 
variations  of  the  derivatives  of  displacements  are  elimi¬ 
nated  through  the  use  of  integration  by  parts,  Eq.  (2.18) 
is  obtained  (written  in  terms  of  stress  resultants). 


sin  a  d6  (2.18) 
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II 

III 

or  u  =  0 

(2.23) 

Nse  -  Ns0  H  T 

or  v  =  0 

(2.24) 

j  f (sMs) *  -  He 

-  2MgQ  esc  a] 

» 1  £s%  ■ 

esc  a  M~g  =  S’ 

1! 

fs 

O 

0  (2.25) 

M_  =  M„  b  M 
s  s 

i 

or  w  =0 

(2.26) 

where  N,  T,  S,  and  M  are  the  generalized  boundary  forces 

i 

associated  with  u,  v,  w,  and  w  ,  respectively.  General¬ 
ized  boundary  forces  here  are  both  forces  and  moments  per 

unit  length  along  the  shell  edge.  The  indicated  deriva¬ 
tives  are  defined  by  the  notation 

in 

/'ol'O 

(2.27) 

If  =  (  > 

(2.28) 

For  the  free-free  case  of  no  displacement  constraints  and 
no  external  edge  forces,  the  boundary  conditions  become 
simply 


N„  =  0  at  s  =  s,  and  s„ 

JL 

(2.29) 

N  .  =  0  at  s  =  s,  and  s0 

G  0  12 

(2.30) 
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[  (sMg ) '  -  Mg  -  2Msg  csca]  =  0  at  s  =  s1  and  (2.31) 

Mc  =  0  at  s  =  and  s0  (2.32) 


For  free  vibrations  the  external  forces  per  unit 
area,  F  ,  Fy,  and  Fw,  are  replaced  by  inertia  forces  using 
D’Alembert's  principle.  Rewriting  the  equations  of  equil¬ 
ibrium  in  terms  of  displacements  such  that  they  physically 
represent  force  per  unit  area  in  the  longitudinal,  cir¬ 
cumferential,  and  radial  directions,  respectively,  Eqs. 
(2.20)  to  (2.22)  become 


d  [tysii 


!  . 


CSC 


(3-v)  v*  » 

a  -  v  g  J — 2  csc  a  +  u 
s 


+  ligl  esc2  a  +  H-  -  »  .  cot  a 


vw 


2  s2 


5  s2 


>2 

+  —5-  cot  a]  -  ph  — =  0 


3t£ 


(2.33) 
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[(yv)  V*  +^csc2a  +Cl^v)|_- 


Ckv)  2_ 

2  2 


.  (1+v)  u  * _ _  (3-v)  u* 

+  7  - csc  a  +  K^Lo~/  ~o 

2S  d 

^2 

-  Ph  =  0 

at2 


csc  a  -  —  cot  a  csc  a] 


(2.34) 


D  csc  a+^-u  +^2"^2  co-fc  cot  a 

-  s  s  s 

Dh2  r..""  .  r,  w'1.. _ 2  _  ,  w” _ 4  „ 

-  22“  Lw  +  2  - g—  csc  a  +  —jr  csc  a 

s  s^ 


,  2w'"  2w'** 2  _  w"  ,  4w“ _ 2  „ 

H — - —  -  — —x —  csc  a  -  —g-  h — — q-  csc  a 


+  w1 ,  .  d2w  _  n 

-VJ  -  ph  — 5-  -  0 

s3  St2 


(2.35) 


where  the  constant  p  is  the  mass  density.  The  boundary 
conditions,  Eqs.  (2.29)  to  (2.32),  written  in  terms  of  dis¬ 
placements  and  units  of  generalized  force  per  unit  length. 


become 


D[u«  +  v(^4°t^.  +  |L  esc  a)]  =  0 


D^^[v'  -  ~  CSC  a]  =  0 


(2.36) 


(2.37) 


•jg—  [w"  +  w csc2  a  -  w“  esc2  a 

s  s 


+  f -4j  =0 
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(2.38) 


[w"  +  v(^tj-  csc2  a  +  —-)]  =  0 


(2.39) 
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CHAPTER  III 


SOLUTION  OF  THE  EQUATIONS 


Form  of  Solutions 


As  described  briefly  in  Chapter  I,  the  emphasis 
during  the  early  phase  of  this  study  was  to  investigate 
the  possibility  (theoretically)  that  pure  vibration  modes 
of  the  form  found  by  Watkins  and  Clary  in  their  experi¬ 
mental  work  are  possible.  In  order  to  permit  solutions 
wherein  a  different  number  of  nodes  can  occur  at  the  two 
cone  edges ,  consider  complete  Fourier  series  solutions 
of  the  form 


00  (1)  00  (2) 

u  =  [  Z  U  (s)  cos  n0  +  Z  U  (s)  sin  n0]  sin  u>t 
.  n=o  n  n=l  n 

(3.1) 


00  (1)  *  (2) 

v  =  [  £  V  (s)  sin  n0  +  Z  V  (s)  cos  n0]  sin  out 
-  n=l  n  n=o  n 

(3.2) 


“  (1)  00  (2) 

w=(Z  W  (s)  cos  n0  +  £  W  (s)  sin  n0]  sin  u)t 
n=o  n  n=l  n 

(3.3) 
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Regarding  the  functions,  (s),  Vn  (s),  etc.,  to  be 

arbitrary,  this  general  form  of  solution  has  the  property 
of  completeness  and  can  represent  any  possible  mode 
exactly.  Consider  now  the  substitution  of  the  solutions, 
Eqs.  (3.1)  to  (3.3) >  into  the  first  of  the  three  appli¬ 
cable  differential  equations,  Eq.  (2.33).  The  result  is 
given  by  Eq.  (3-4)  . 


.  .  oo  (l)*  05  (  2  )  ’ 

D[(l+¥.)  c.s.ca.  [  e  n  V  cos  n9  +  E  n  V  sin  n0] 

d  s  - n=l  n  n=o  n 

00 

(3-v)  csca  r  ”  ..  ,r  ^  cos  n0  4-  E  „  -r  ^2'  sin  n8] 
'  2  “2“  L  V1  Vn  r;=on  "n 


s  n=l 


+  C  2  U 
-  n=o 


(1)  ’ 1 


cos  n0  +  E  U 


(2)'« 


sin  n0] 


o  00  /  1  }  00  f2) 

.  (!z^)  a  [  2  n2U  cos  n0  +  E  n2  U  sin  n0] 
2  sd  n=o  n  n=l  n 

00  (1)  00  (2) 
f  —  [  2  U  cos  n0  +  £  U  sin  n0j 
3  n=o  n  n=l  n 

1  -  (1)  *  (2) 

-  -=rr  [  E  U  cos  n9  +  E  U  sin  n9] 
s  n=o  n=l 


v  cot  a 


[  E  Wr 


cos  n9  +  E  W„ 


sin  n0] 


cot  a  r  "  ^  “  (2) 

+  .£oua  [  e  w  cos  n0  +  E  W  sin  n9]} 

s^  n=o  n  n=l  n 

2  <0  C1)  *  (2)  ,  ,  . 

+  phuj  [  E  Un  cos  n9  +  E  Un  sin  n9]  =  0  (3*^) 

n=o  n=l 


Equation  (3-4),  plus  two  other  similar  equations 
which  are  obtained  when  the  displacements,  Eqs.  (3-1)  to 
(3-3) )  are  substituted  into  the  remaining  equations  of 
motion,  Eqs.  (2.3^)  and  (2.35) >  are  to  be  solved  for  the 

U)  ,  W,  s 

eigenfunctions  Un  (s'>  (s)>  etc.  and  the  cor- 

2 

responding  eigenvalues  uj  .  Following  the  usual  procedure, 
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these  three  equations  are  multiplied  by  either  cos  m9  or 
sin  m0  and  integrated  from©  =  0  to  0  =  2tt.  Consider,  for 
example,  the  multiplication  of  Eq.  (3-^)  by  cos  m0  and  the 
subsequent  integration  from  0  =  0  to  0  =  2ir.  Due  to 
orthogonality,  the  terms  involving  the  functions  with  a 
superscript  (2)  are  completely  eliminated,  the  only 
remaining  terms  being  those  for  which  n  =  m  and  involving 
the  superscript  (1) .  When  both  integrations  are  performed 
independently  on  each  of  the  three  equations ,  the  result  is 
six  equations,  three  of  which  involve  the  functions 
U  (s ) ,  vm^(s)>  and  Wm^1^  (s)  while  the  remaining  three 
involve  the  functions  Ura^(s)>  and  Wm(2^(s). 

Furthermore,  the  three  equations  involving  the  functions 
with  superscript  (1)  are  identical  with  the  three  equations 
involving  the  superscript  (2).  Thus  there  is  no  need  to 
consider  both  sets  of  functions.  The  family  of  modes  for 
which  n  =  m  is  completely  given  by 


u  =  U  (s)  cos  m9  sin  out 
mv  ' 


v  =  V  (s)  sin  m0  sin  out 
mv  ' 


w  =  W  fs)  cos  m9  sin  out 
m'  ' 


(3.5) 

(3.6) 

(3.7) 


where  m  is  the  number  of  full  waves  around  the  circum- 


ference  of  the  shell,  also  defined  by  one-half  the  number 
of  nodal  points  encountered  when  going  around  the  circum¬ 
ference. 

For  the  purpose  of  obtaining  an  approximate  solu¬ 
tion  through  the  use  of  the  Galerkin  procedure,  a  finite 
series  approximation  of  the  form 

n 

u  =  £  q.  a.(s)  cos  m0  sin  rnt 
i=l  1  1 

n 

v  =  £  q.  p.(s)  sin  m9  sin  uut 
i=l  1  1 

n 

w  =  £  q.  Y-(s)  cos  m0  sin  tut 
i=l  1  1 

was  used  for  this  investigation.  Here  m  is  allowed  to  take 
on  a  single  value  for  a  specific  solution.  The  choice  of 
the  functions  ai(s),  P±(s)  and  (s )  are  discussed  further 
in  Chapter  IV.  Modal  patterns  resulting  from  using  these 
displacements  will,  of  course,  consist  of  parallel 
circles  and  equispaced  meridians  —  patterns  which  because 
of  the  axial  symmetry  of  the  shell  and  the  boundary  condi¬ 
tions,  and  the  circumferential  periodicity  of  the  vibra¬ 
tion  motion,  are  entirely  logical.  The  experimental 
investigation  of  Hu,  Gormley,  and  Lindholm  (6),  conducted 
at  the  Southwest  Research  Institute,  confirms  the  foregoing 


(3.8) 

(3-9) 

(3.10) 
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discussion.  Further  details  of  these  results  appear  in 
Chapter  V. 


Application  of  the  Galerkin  Procedure 

The  Galerkin  procedure  used  by  Bolotin  (1)  and 
Matthews  (8)  differs  from  the  usual  Galerkin  procedure  in 
that  the  functions  utilized  for  the  assumed  solution  must 
satisfy  only  the  displacement  or  "natural"  boundary  condi¬ 
tions.  This  is  an  important  feature  since  in  the  case  of 
the  free-free  conical  shell  it  is  inconvenient  to  choose 
functions  which  satisfy  the  boundary  conditions,  Eqs. 
(2.36)  to  (2.39).  In  fact,  it  is  difficult  to  find  func¬ 
tions  which  do  satisfy  these  boundary  conditions.  In  this 
method  the  errors  in  the  forces  and  moments  at  the  bound¬ 
aries  are  weighted  along  with  the  error  functions  from  the 
differential  equations.  The  total  weighted  error  is  then 
set  equal  to  zero  which  is  equivalent  to  setting  the 
virtual  work  done  by  all  the  forces  and  moments  moving 
through  appropriate  displacements  equal  to  zero. 

The  modified  Galerkin  procedure  requires  addi¬ 
tional  equations  to  account  for  the  boundary  errors.  The 
equations  are  obtained  from  the  requirement  that  boundary 
forces  and  moments  must  balance  the  internal  forces  and 
moments  on  an  infinitesimal  end  element.  Utilizing  the 
force  type  boundary  conditions,  Eqs.  (2.36)  to  (2.39),  in 


where  N,  T,  S,  and  M  are  the  generalized  boundary  forces 
per  unit  length  as  previously  defined.  Similarly,  the  set 
of  equations  for  the  other  edge  of  the  shell  where  s  =  s^ 
is  given  by  Eqs.  (3.15)  to  (3.18).  In  Eqs.  (3*11)  to 
(3.14)  above  as  well  as  in  Eqs.  (3-15)  to  (3.18)  below, 
positive  forces  are  generally  taken  to  be  in  the  direction 


23 

of  positive  displacements  (see  Figures  2.1  and  2.2).  Care 
must  be  exercised  in  the  determination  of  the  signs  of  the 
forces  at  the  shell  boundaries.  The  signs  which  appear  in 
these  equations  must  be  carried  through  the  steps  of  the 
modified  Galerkin  procedure  which  follow. 


! 

I 


(N  -  D  [u1  +  g-  (u  -  w  cot  a)  4-  — esc  a]} 

s 


=  o  (3.15) 

3=  S^ 


{T  -  [v1  -  esc  cij} 

s  s 


=  o 
s=  sr 


(3*16) 


(s  +  ^-  [w'"  +  g )  w'“  csc2a  -  W  csc2a 

s^  S'5 


"  t 

,  w  W1-,-) 

+  ~‘71! 


=  0 
s=  s. 


(3-17) 


{-  M  -  [w"  +  v(^g-  csc^a  +  “-*) ] } 


2  .  w 1 


r° 

s=  s. 


(3.18) 


An  exact  solution  must  satisfy  the  end  element 
equations,  Eqs.  (3.11)  to  (3.18),  as  well  as  the  general 
element  differential  equations,  Eqs.  (2.33)  to  (2.33),  The 
approximate  solution  to  be  considered  here  will  not  satisfy 


I 


Jr*#-  r*.  -« 
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all  of  these  equations.  Therefore,  errors  in  all  of  them 
are  to  be  accounted  for  in  the  modified  Ga.lerkin  procedure. 

A  set  of  error  functions  is  generated  by  the  sub¬ 
stitution  of  the  assumed  displacements,  Eqs.  (3.8)  to  (3*10), 
into  the  general  equilibrium  equations,  Eqs.  (2.33)  to  (2.35) 
and  the  boundary  equilibrium  equations,  Eqs.  (3.11)  to  (3.18). 
These  are  given  by  Eqs.  (3*19)  to  (3*29)  which  follow. 


cos  m0  sin  out  E  {D[ 

(3-v)  esc  a  „  ( 1— v )  esc  a  „2„  ,  1  „  , 

-  Wg-'  - +  a±  -  '-g-'  —5—  m  ai  +  T  ai 


1  _  v  cot  a  ,  ,  .  cot  a 

'  JL  CC  .  —  ■  ■  -  17  1  -1-  - 

7?  i  s 


v1'  +  COt£g  Yi J  +  phu-'2ai) 


=  «a(s.8) 


2 

sin  m0  sin  cot  E  q±  {D[^p)  P±'  -  CSC'2a  m2p1 
i  -  s 

+  (l=yj  8  1  -  (^v)  e  .  (1+v )  csc_a  t 
+  2s  Pi  2s2  pi  2  s  m  ai 

(3-v)  esc  a _ _  cot  a  esc  a  _  .. 

~  Wr*'  - 2~  in  a.  + - ^ - -  m  YiJ 


(3.19) 


+  phui20i}  =  €p(s,0) 


(3.20) 


~\*r=  --  - 
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cos  m9  sin  tut  E  q.  {D[-s--^-  mpi  +  ^  ai* 
i  s 


v 


p 

.  1  „  cot  a  .  ,  ,  Dh  r  >“' 

+  —  a± - Yi3  cot  a  -  -jg-  [y± 

S  S  -  a 

+  ^  >4v±  +  |  Vi'  4-  SSgA  m2Yi'  -  %  Vi' 


2csc  a  _2..  « 
- *—  m  y± 


4cscda  2.  .  1 


Tf—  m  Y^  +  —j  Y^'J  +phu)Y^}  =eY(s,0)  (3.21) 


-  N  +  cos  m8  sin  tut  E  q.  {D[a.  ’  +  —  (a.  -  y.  cot  a 

^  1  X  SI  1 


+  esc  a  mPi)]3 


■  Ei<9> 


(3.22) 


3=S. 


N  -  cos  m0  sin  tut  S  q.  {D[a.  *  +  ~  (a .  -  y.  cot  a 

^  X  X  S  X  X 


+  esc  a  mg^)]} 


-  V0) 


(3.23) 


S=S. 


it-  « 


•,~SS>i=rf 


r 


^«steSKs?»i4 


26 


-  H  +  sin  m0  sin  cut  £  qi  {.DC^“.V)  [p^i  -  ~ 


esc  a 
s 


mai]} 


=  ^(e) 


(3.24) 


s=s. 


T  -  sin  m0  sin  tut  £  q^  tP^'  “  ir  Pj[ 


esc  a 


ma±] } 


E4(6) 


(3.25) 


s=s. 


2 

-  S  -  cos  me  sin  cut  E  q±  [Yj*  -  csc2a  m2Yi’ 


+  Ql^.)  Csc2a  m2Y1  +  — | - -^r  Y^  *  ] } 


=  E5(0)  (3-26) 


s=s. 


2 

■§■  +  cos  me  sin  art  E  q±  [ Y±w  -  ^5^  csc2a  m2Yi’ 


+  csc2a  m2Yi  +  — | ^2  Y± *  ] } 


=  Eg(6)  (3.27) 


s=s. 


1 


M 
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M  +  cos  m6  sin  u>t  2 


esc  a  _2 


r,  _  run  r..  »  ...  /  esc  i 

5  qi  lI5“tYi  +  v  ( - 75 

1  s 

iYi')))  =^(8) 


m  v. 


(3.28) 


's=s. 


2  2 

H  -  cos  ni8  sin  lot  2  q^  [y^  +  v  (-  c-s-?~  m2vi 

i  s 


tV»> 


-  Eg(e) 


(3.29) 


s=ss 


In  the  modified  Galerkin  procedure,  the  error 
functions  representing  the  end  errors,  Eqs.  (3.22)  to 
(3*29) ,  are  weighted  by  displacements  in  a  manner  similar 
to  the  weighting  of  the  error  functions,  Eqs.  (3*19)  to 
(3.21),  resulting  from  the  differential  equations  of 
motion.  Error  functions,  ea(s,8),  E^(8),  and  Eg(8),  which 
represent  errors  in  force  in  the  longitudinal  direction, 
will  be  multiplied  by  the  displacement  u  per  unit  dis- 
placement  of  the  «j  coordinate.  Similarly,  e^(s,8), 
E^(8),  and  E^(8),  which  represent  errors  in  force  in  the 
circumferential  direction,  will  be  multiplied  by  the  dis- 
placement  v  per  unit  displacement  of  the  j  coordinate; 


r 


m,  .»  M  jmij  iwa  ' '  ii.il.il 


*^*-^*»-  V  >C*ac&  -  #v>rs.  v^= 
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and  €y(s,0),  E^(6),  and  Eg(@),  which  represent  errors  in 

force  in  the  radial  direction,  will  he  multiplied  by  the 

th 

displacement  w  per  unit  displacement  of  the  j  coordinate. 
Functions  E^(0)  and  Eg(0),  which  represent  errors  in 
moment,  will  be  multiplied  by  the  rate  of  change  of  radial 
displacement  in  the  longitudinal  direction  dw/ds  per  unit 
in  the  coordinate.  The  sum  of  all  these  weighted 
errors,  integrated  over  the  er.tire  conical  shell,  will  be 
set  equal  to  zero.  Written  in  terms  of  the  error  func¬ 
tions  defined  previously,  application  of  the  foregoing 
procedure  yields  the  equation,  Eq.  (3-30) • 


27 r  s. 


J  J  l-eaaj  cos  m0  +  eg8j  Sin  m6  +  eyYj  COS  m0^S  Sin  a 


dGds 


o  s. 


27T 


+  J  fExa<j  cos  3318  j  +  E2aj  cos  m8 

-  S=S  -i  S  =S  f 


+  j  sin  m0  |  +  E^p ^  sin  m8 

S^Sj 


+  Ej-Yj  cos  m0 
s=s2  5  J  's=s1 


+  E/y.  cos  m9 |  +  E  y,1  cos  m0 

6  J  s=Sg  7  J 


s=s. 


+  EgY j 1  cos  m9  I  ] s  sin  a  d9  =  0 

S=S  r\ 


(3.30) 
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The  foregoing  procedure  is  equivalent  to  applying  the 
principle  of  virtual  work,  i.e.,  requiring  that  the  vir¬ 
tual  work  done  by  the  forces  and  moments  in  a  displacement 

•f* 

in  the  j  coordinate  be  equal  to  zero.  Expanding  Eq. 
(3-30)  through  the  use  of  the  error  functions,  Eqs.  (3.19) 
to  (3.29),  integrating  the  result  over  8  and  cancelling 

i.u  xu 

terms,  the  general  i  n  term  in  the  J  equation  can  be 
expressed  as 


J  2  ((l+v)  S2SJL  ^  .  (Sgv)  =50^  mpl0j  + 

S1 

p 

(1-v)  CSC  a  _2_  _  .  1  _  V  cot  a  ..  1  _  _ 

—  ' — 75—*^  '  '  '  ®  j  +  ~  Cc^ 1  Oc  j  —  ^  y ^ '  cij  —  — ^  u _^d . 


.  cot  a  _  .  ^(l-v^)  ...2 _ _  (1-v)  „  » n 

+  — +  p  — w  a±aj  +  —  P±  P 3 


CSC  a  2  (1-v)  o  ,R  (1-v)  o  a 

s- m  H*S  +  <nr  pi  ej  -  'jp  eipo 


(1+v)  esc  a _ tr>  (3-v)  esc  a _ _ 

-  — r-”ai ■»}  - 


cot  a  esc 


-  mYjPj  +  "SjPj  +  mpiY- 


V  t  1  Cot)  (X  1  ,  .  ».  w/ 

+  8  °i  YJ  +  7?  aivj  -  —5-  YiYj]  1  -  k  tY±  Yj 

s  s 
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2  2i 

2  esc  a  _2_.  ///..  .  esc  a  _4 . 2  .. 

- <5 —  m  Yi  +  — n?-  m  YiYj  +  s'  Yi  Yj 


,  2  esc2  J2...  1  ..  //..  4  csc2a  _2..  .. 

+ - ^ —  m  Yi  Yj  "  *"5"  Yi  Yj  "  - 4 —  m  YiYj 

s  s  s 


2 

+ ' Yj 3  +  ^  sin  a 


ds 


+  [a^Oj  +  vaia«1  - 


v  cot  a  v  esc  a  „  i 

- - - "  a.  + - - —  mp.a.Js  sm  a 

S  J.  J  s  1  J 


-  -  -H  -  Yiaj  -  2-£Sra  ^  *  I 


+  tpi'pj  "  J  "  £fifJL  na^jls  sin  a 


-  )  rptp  Ipp  esc  a 

2^  lPi  Pj  7  pipj - I 


ma^Pj] s  sin  a 


+  ktYi'V  ~  ~  °°l  a  m2YiYj '  +  T  Yi’Vls  sln  a  | 


*  [v/Yj1 


p 

v  esc  a  _2  i  .  v  ,  , t 

- 2 —  m  YiYj  +  J  Yi  Y j  Is  sin  a 
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-  Qgii)  s  p1»pJ'  +  [p^^P  »28  -  52SJL 


jggJ^JL  (m2csc2a  -  4)]  y^Yj  -  -|(2m2csc2a  +1)  Yi  *  v ^  * 


-  ksVY/  -  ^  («*  +  Vj) 

+  m  CSC  a  (a±pj '  +  P^a^)  +  (a±y^  +  y^j) 


,  /  i  f  \  #  m  CO  u  GC  CSC  CC  /  0  ,  0  \  *1 

+  V  cot  a  (ai’y.  +  y^’)  + - g -  (y±P j  +  P^j)} 


ds 


.  r  ...  -  .  ( 1-v)  „  „  .  km2(3-v)  csc2a  ..  .. 

+  t-  va^,  +  '-5-'  H -  >iY; 


+  kvm2cscgg  (Y±Yj.  +Yl>Yj)  -  fcvYi'Yj' 


-  vm  esc  a  (a.^^  +  P.^)} 


(3-33) 


Equation  (3-33) >  which  is  the  basis  for  all  calculations 
that  follow,  may  also  be  derived  using  fche  Rayleigh-Ritz 
procedure  beginning  with  an  energy  expression  of  the  form 
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C-U  *2 

o  -  J  J'  I  CD  (€s2  +  €e2  +  2VSs€e  + 

O  S.^ 

+  kD  (Ks2  +  Ke2  +  2vKgKQ  +  2(l-v)  Kg02 
-  (Fuu  +  Fyv  +  Fww)]  s  sin  a  dsd0 


(1= v) 
2 


(3.34) 


and  utilizing  the  strains,  Eqs.  (2.1)  to  (2.3) >  and  the 
curvatures,  Eqs.  (2.7)  to  (2.9).  Details  of  this  deriva¬ 
tion  will  not  be  presented  here.  The  Galerkin  and 
Rayleigh- Ritz  procedures  are  equivalent  when  the  differen¬ 
tial  equations  of  equilibrium  used  in  the  Galerkin  method 
are  obtained  from  a  variation  of  the  energy  expression 
used  in  the  Rayleigh- Ritz  procedure.  This  was  pointed  out 
by  Singer  (11). 


CHAPTER  IV 


COMPUTATIONAL  PROCEDURE 
Choice  of  Functions 

For  the  case  of  a  conical  shell  free  of  edge  sup¬ 
ports,  the  choice  of  displacement  functions  in  the  longi¬ 
tudinal  direction,  ai(s),  g^(s)  and  Y^(s),  are  somewhat 
arbitrary  since  the  displacement  boundary  conditions  are 
completely  relaxed.  Polynomials  have  been  used  by  sev¬ 
eral  investigators,  i.e.,  Platus  (9),  for  the  case  of 
fixed-end  conditions  on  one  cone  edge.  Because  of  the 
great  number  of  integrations  that  must  be  performed  to 
determine  the  matrix  coefficients,  it  was  decided  to  also 
utilize  polynomials.  There  are,  of  course,  other  possi¬ 
bilities  such  as  the  functions  which  might  be  obtained  from 
the  free  vibrations  of  a  tapered  beam  without  end  supports. 
However,  the  necessary  integrations  would  be  unwieldy. 

Third  degree  polynomials  of  the  form  of  Eqs. 

(4.1)  to  (4.3)  were  assumed  in  the  determination  of  the 
results  presented  in  this  study.  Higher  degree  polyno¬ 
mials  are  better  able  to  represent  the  more  complex  shapes 
of  modes  higher  than  the  third.  On  the  other  hand,  the 
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size  of  the  matrix  which  must  be  inverted  is  increased  by 
three  rows  and  columns  for  each  added  degree,  a  factor 
which  reduces  computational  accuracy.  Primary  interest  was 
in  the  first  two  modes  for  each  case  of  a  fixed  number  of 
circumferential  waves.  Therefore  a  degree  high  enough  to 
adequately  represent  the  lower  modes  was  determined  to  be 
the  logical  choice  for  this  investigation.  Further  dis¬ 
cussion  on  computational  accuracy  is  given  later  in  this 
chapter. 

o  i-1 

a.  =  (•§—)  i  =  1  to  4  (4.1) 

1  S2 

=  (f-)1  5  i  =  5  to  8  (4.2) 

1  S2 

.  i-9 

Yj  =  (|— )  i  =  9  to  12  (4.3) 

i  s2 

Other  functions,  i.e.,  a^,  a,g,  p,,  pg,  etc.,  were  taken  to 
be  zero.  This  seemingly  odd  indexing  was  necessary  in 
order  to  assemble  the  matrix  elements. 


Arranging  the  Equations  for  Computer  Solution 


The  n  homogeneous  equations  which  result  from  the 
application  of  the  modified  Galerkin  procedure  to  the 
problem  of  free  vibrations  of  a  thin  truncated  conical 
shell  can  be  expressed  in  the  matrix  form 


P 


(i=£) 

E 


u)2[M] 


pi’ 

1 

[>l 

. 

[V 

'  -  [K]  - 

1  is 
'  <■ 

=  • 

(4.4) 


where  the  generalized  mass  matrix  [M] ,  determined  from  the 
general  ijth  term,  Eq.  (3-33)>  is 

s2 

Mfj  =  J  (aiaj  +  PjPj  +  YiYj)  sds  (^*5) 

S1 

and  the  generalized  stiffness  matrix  [K] ,  also  determined 
from  Eq.  (3 . 33 ) ,  is 


- 
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Having  chosen  the  functions  a^,  and  y^,  the  above  mass 
and  stiffness  matrix  elements  can  be  evaluated  for  a  shell 
having  specified  geometry  and  a  fixed  number  of  circumfer¬ 
ential  waves,  m. 

In  order  to  solve  t...  problem  numerically  on  a 
ligital  computer,  matrix  Eq.  (4.4)  can  be  more  favorably 
rewritten  as 


p  (1-v2)uj2/E 


r  ~  ■'i 

ql 

• 

• 

• 

►  =  [K]"1  EM]  , 

>1 

K. 

1 

[x 

(4-7) 


where  [K]  1  is  the  inverse  of  the  matrix  [K] .  Equation 
(4.7)  is  solved  through  a  process  of  matrix  iteration 
which  produces  the  n  frequencies  and  corresponding  mode 
shapes  for  the  simplified  constrained  system.  After  the 
fundamental  frequency  and  mode  shape  is  found,  a  process 
of  matrix  sweeping  is  utilized  to  remove  the  first  mode, 
thereby  allowing  the  second  iteration  process  to  converge 
on  the  second  mode.  This  technique  is  continued  until  as 
many  frequencies  and  mode  shapes  as  are  desired  have  been 
determined,  up  to  a  maximum  of  n.  Computations  were 
accomplished  on  the  IBM  7072  computer  at  The  University  of 
Arizona  and  the  IBM  7094  computer  at  the  Western  Data 
Processing  Center,  University  of  California  at  Los  Angeles. 
The  program  is  similar  to  that  used  by  Matthews  (8) 
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although  somewhat  more  complex.  The  FORTRAN  language  pro¬ 
gram  consisting  of  the  main  program  and  subroutines  for 
matrix  inversion,  matrix  multiplication  and  matrix  itera¬ 
tion  are  presented  in  Appendix  B. 

Accuracy  of  Computational  Procedure 

The  formulation  of  the  equations  is  such  that  the 
matrix  which  must  be  inverted  approaches  singularity  as  the 
half  cone  angle  approaches  zero.  This  is  due  to  the  fact 
that  the  equations  of  equilibrium,  Eqs.  (2.33)  to  (2.35), 
degenerate  as  a  approaches  zero  and  a  limiting  process 
must  be  utilized  for  the  case  a  =  0.  In  addition,  for  any 
reasonable  length  cone,  the  quantities  s-^  and  s2  become 
very  large  as  the  cone  angle  is  made  small.  Since  differ¬ 
ences  of  high  powers  of  s^  and  s2  are  needed  in  the  compu¬ 
tation  of  the  matrix  coefficients,  round-off  errors  and 
loss  of  significant  figures  degrade  the  results.  As  a  con¬ 
sequence  of  these  factors,  computed  frequencies  and  mode 
shapes  for  small  cone  angles  are  expected  to  be  less 
accurate  than  those  for  larger  cone  angles. 

Terms  in  the  expressions  for  curvature,  Eqs.  (2.4) 
to  (2.6),  which  are  neglected  in  Eqs.  (2.7)  to  (2.9)  become 
less  significant  as  a  approaches  90°  since  the  terms  contain 
cos  a.  This  fact  also  contributes  to  greater  accuracy  for 
the  larger  cone  angles. 
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Second  degree  polynomials  of  the  general  form  of 
Eqs.  (4.1)  to  (4.3)  were  utilized  for  certain  cases  to 
determine  the  effect  of  choice  of  degree  of  the  assumed 
longitudinal  displacement  polynomials  upon  the  results. 
Differences  betw  n  computed  frequency  parameters  were 
found  to  increase  as  the  circumferential  wave  number  m  is 
decreased  and  as  the  mode  number  is  increased.  For 
example,  in  a  typical  case,  the  difference  in  frequencies 
for  the  first  mode  is  about  two  percent  for  m=2  and 
decreases  to  less  than  one-tenth  of  a  percent  at  m=9> 
whereas,  for  the  third  mode  the  difference  for  m=2  is  nine 
percent  and  decreases  to  one  percent  at  m=9*  Table  I 
illustrates  this  point,  wherein  the  results  are  normalized 
such  that  the  frequencies  determined  utilizing  the  four 
term,  third  degree  polynomials  are  equal  to  one,  thereby 
allowing  convenient  comparison  with  frequencies  determined 
from  the  use  of  the  second  degree  polynomials. 

As  m  increases,  the  complexity  of  the  circumfer¬ 
ential  shape  appears  to  compensate  for  the  lack  of 
accuracy  in  the  longitudinal  shape  of  the  mode.  That  is, 
for  the  higher  circumferential  wave  number  shapes,  the 
frequency  is  predominately  determined  by  the  complexity 
of  the  circumferential  shape.  The  inaccuracy  of  the  lon¬ 
gitudinal  shape  does  not  seriously  effect  the  natural 
frequencies  of  the  lower  modes.  On  the  other  hand,  the 


TABLE  I 

NORMALIZED  FREQUENCIES  FOR  THREE  AND  FOUR  TERM  POLYNOMIALS 


1st  Mode 

2nd  Mode 

3rd  Mode 

m=2 

Three  Term 

102.2 

105.4 

108.9 

Four  Term 

m=4 

Three  Term 

Four  Term 


100.0 


100.3 

100.0 


m=9 

Three  Term 

Four  Term 


100.03 

100.00 


100.0 


105.3 

100.0 


100.05 

100.00 


100.0 


108.8 

100.0 


101.0 

100.0 


inability  of  the  less  complex  (lower  degree  polynomial)  lon¬ 
gitudinal  assumed  shape  to  adequately  represent  the  more 
complex  modes  is  evident,  i.e.,  the  higher  degree  polyno¬ 
mial  is  better  able  to  represent  the  higher  modes,  which 
is  to  be  expected. 
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In  this  study,  the  prd-mary  interest  is  in  the  lower 
modes,  especially  the  first  and  second.  Therefore,  it  was 
decided  to  utilize  the  third  degree  polynomials  as  given 
by  Eqs.  (4.1)  to  (4.3)  for  the  longitudinally  varying  dis¬ 
placement  function.  The  use  of  higher  degree  polynomials 
would  not  significantly  increase  the  accuracy  of  the  fre¬ 
quencies  and  mode  shapes  of  the  first  two  modes,  based 
upon  results  for  the  test  case  using  both  second  and  third 
degree  polynomials. 

According  to  Bolotin  (1),  convergence  in  connec¬ 
tion  with  the  Galerkin  procedure  has  been  investigated 
only  in  connection  with  very  simple  problems.  The  wide  use 
of  Galerkin1 s  method  for  problems  in  which  convergence  is 
not  conclusively  proved  is  based  upon  a  comparison  of 
existing  exact  results  with  those  obtained  with  Galerkin’ s 
method,  for  those  few  cases  where  exact  solutions  can  be 
obtained  analytically.  One  must  rely  to  a  certain  extent 
on  checks  built  into  a  computational  procedure  which  may 
be  designed  to  identify  non- convergence.  A  check  on  mode 
orthogonality  with  respect  to  the  generalized  mass  matrix 
was  utilized  in  this  study. 
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CHAPTER  V 

COMPARISON  OF  ANALYTICAL  RESULTS  WITH  THE  EXPERIMENTAL 
DATA  OF  HU,  GORMLEY,  AND  LINDHOLM 

Geometry  of  the  truncated  conical  shells  taken 
from  Ref.  6,  and  used  here  as  comparative  cases,  is  given 
in  Table  II.  Parameters  are  as  shown  diagrammatic ally  in 
Figure  2.1. 

TABLE  II 

CONE  GEOMETRY,  CASES  CF  HU,  GORMLEY,  AND  LINDHOLM 


Case  1 

Case  2 

Case  3 

Case  4 

Cone  Half- 
Angle,  a 

14.24° 

30.24° 

45.12° 

60.42° 

Major  Edge 
Radius,  rL 

6.07  in. 

7.95  in. 

8.96  in. 

10.00  in. 

Minor  Edge 
Radius ,  rQ 

2.72  in. 

3.49  in. 

3.98  in. 

4.45  in. 

Length,  L 

13.2  in. 

7.65  in. 

4.96  in. 

3.15  in. 

Thickness,  h 

0.01  in. 

0.01  in. 

0.01  in. 

0.01  in. 

Completeness 

Parameter, 

Sg/Si 

2.23 

*  2.28 

2.24 

2.25 

43 


Comparative  results  for  Cases  1,  2,  5  and  4  are 
given  by  Tables  III  through  VI.  The  frequency  parameter 
0  for  these  cases  is  defined  by 

(1  2,  1/2 

n  =  u,rL[P  E  ]  (5-1) 

Agreement  between  the  theoretical  and  experimental  data 
appears  to  improve  as  the  cone  angle  increases.  This 
trend  is  partially  due  to  the  fact  that  computational 
errors  increase  as  the  cone  angle  approaches  zero,  as 
explained  in  Chapter  IV.  In  general,  the  agreement  also 
improves  in  each  case  as  m  increases.  Due  to  the  low 
frequencies  being  measured,  the  experimental  results,  of 
course,  were  air  o  subject  to  error.  Errors  introduced  by 
cone  imperfections,  seams,  supports,  and  method  of  excita¬ 
tion  were  most  likely  greater  for  the  lower  circumferen¬ 
tial  wave  numbers. 

It  can  be  seen  that  the  computer  results  agree 
quite  well  with  the  experimental  data.  It  can  also  be 
noted  from  Tables  III  to  VI  that  with  few  exceptions  the 
theoretical  frequencies  are  grea.ter  than  the  experimentally 
determined  ones .  Again  this  is  to  be  expected  because  of 
the  constraints  imposed  by  an  approximate  method. 

Mode  shapes  for  Case  3  (45-12°  cone  half  angle) 
are  presented  in  Figures  5.1  and  5-2.  Actually  the  nor¬ 
malized  radial  displacement  is  plotted  rather  than  the 
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TABLE  III 

FREQUENCY  PARAMETERS,  CASE  1 


Circum- 

Frequency  Parameter  0 

XcicI/OiaJ.  ” 

Wave 

1st 

Mode 

2nd 

Mode 

No. 

NASA- CR- 384 

Galerkin 

NASA- CR- 384 

Galerkin 

2 

.00191 

.00243 

.0113 

.0062 

3 

.00502 

.00571 

.0195 

.0144 

5 

.0151 

.01585 

.0393 

.0394 

7 

.0285 

.0300 

.0611 

.0650 

12 

.0763 

.0789 

.112 

.1285 

18 

.165 

.169 

- 

•=t 
v 0 

OJ 

• 

TABLE  IV 

FREQUENCY  PARAMETERS 

,  CASE  2 

Circum¬ 
ferential  - 
Wave 

Frequency  Parameter  Q 

1st  Mode 

2nd 

Mode 

No. 

NASA- CR- 384 

Galerkin 

NASA- CR- 384“ 

Galerkin 

2 

.00151 

.00191 

- 

.00578 

3 

.00422 

.00456 

- 

.0128 

5 

.0122 

.0126 

- 

.0333 

7 

.0242 

.0244 

.0626 

.0627 

12 

.0647 

.0656 

.120 

.129 

18 

•  135 

.1356 

.200 

.213 
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TABLE  V 

FREQUENCY  PARAMETERS,  CASE  3 


Circum- 

Frequency  Parameter 

n 

X  G»  X*  0  XI  v  1  o<  X 

Wave 

1st 

Mode 

2nd 

Mode 

No. 

NASA- CR- 384 

Galerkin 

NASA- CR- 384 

Galerkin 

2 

.00147 

.00169 

- 

.00612 

3 

.00399 

.00413 

- 

.0130 

5 

.0113 

.01143 

- 

.0313 

7 

.0218 

.0221 

.0583 

.0578 

12 

.0399 

.0607 

.127 

.133 

18 

.123 

.124 

.187 

.212 

TABLE  VI 

FREQUENCY  PARAMETERS,  CASE  4 


Circum¬ 
ferential  - 
Wave 

No. 

Frequency  Parameter 

n 

Is 

t  Mode 

2nd  Mode 

NASA- CR- 3 84 

Galerkin 

NASA-CR-384 

Galerkin 

2 

.00139 

.00144 

- 

.00663 

3 

.00377 

.00367 

- 

.0133 

5 

.0104 

.0104 

- 

.0298 

7 

.0200 

.0200 

.0573 

.0533 

12 

.0552 

.0548 

•  123 

.1245 

18 

.113 

.1125 

.187 

.196 

total  displacement ,  Since  the  radial  displacement  is  dom¬ 
inant  for  the  first  two  modes,  it  was  cnosen  for  discus¬ 
sion  purposes.  Hereafter  the  expression  "mode  "hape”  is 
taken  to  mean,  the  normalized  radial  displacement  shape  when 
reference  is  made  to  the  figures .  In  general,  the  results 
agree  favorably  with  those  found  experimentally  by  Hu  and 
presented  in  Ref.  6.  Figure  5.3  shows  a  comparison  of  the 
experimental  and  theoretical  mode  shapes  for  the  first  and 
second  modes  of  Case  3  for  m=12. 

In  Ref.  6,  Hu,  et  al. ,  propose  the  following  semi- 
empirical  frequency  equation  for  thin  free-free  conical 
shells : 

n  - - ^  m  fa-1)  (rn  +  1  +  4  sin  -S2)  (5.e) 

“n*  -/7TT  2 

It  was  shown  that  agreement  between  the  experimentally 
determined  frequency  parameter  and  that  given  by  Eq.  (5*2) 
was  in  general  very  good  for  the  range  of  parameters 
tested.  It  was  also  stated  that  the  range  of  applicability 
of  Eq.  (5*2)  with  regard  to  the  completeness  parameter 
s2/si  is  uncertain.  Equation  (5.2)  was  developed  to  com¬ 
pare  with  experimental  results  having  a  completeness 
parameter  of  approximately  2.25. 


Figure  5-3  Comparison  of  Shapes 
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In  order  to  determine  the  effect  of  completeness 
parameter  upon  frequency  parameter,  the  basic  shell  of  Case 
3  was  utilized,  varying  sg  while  holding  constant. 
Results  of  varying  Sg/s^  from  1.1  to  10  are  presented  in 
Figure  3-^.  It  can  be  seen  that  for  m=6  and  above,  the 
total  variation  in  frequency  parameter,  for  the  range  of 
Sg/s ^  considered,  is  almost  constant.  This  means  that  the 
effect  of  completeness  parameter  becomes  less  pronounced 
as  m  increases,  as  one  might  suspect.  For  example,  the 
approximate  percentage  increace  in  frequency  parameter 
from  Sg/s^  ~  1*1  to  Q^sl  =  *or  Is  ^5  percent 
while  at  m=12  it  is  only  20  percent.  It  should  also  be 
noted  that  the  variation  in  frequency  parameter  between 
Sg/s-  =1.5  and  Sg/s-^  =  10  amounts  to  about  half  of  the 
total  variation  for  the  range  of  Sg/s^  considered.  Since 
the  Sg/s^  =  2.25  curve  very  nearly  divides  the  area 
between  the  Sg/s^  =  1.5  and  Sg/s^  -30  curves,  reasonably 
good  results  can  be  expected  of  Eq.  (5.2)  for  very  thin 
shells  having  a  completeness  parameter  between  1.5  and  10, 
especially  for  the  larger  values  of  m. 

Attempts  were  made  to  incorporate  the  completeness 
parameter  Sg/s^  Into  Eq.  (5-2),  thereby  extending  its 
range  of  applicability .  However,  these  attempts  were 
unsuccessful.  A  complete  parametric  study  on  complete¬ 
ness  parameter  and  Its  effect  in  relation  to  cone  angle. 
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Figure  5*^  Effect  of  Varying  Completeness  Parameter 


CHAPTER  VI 


COMPARISON  OF  ANALYTICAL  RESULTS  WITH  THE  EXPERIMENTAL 
DATA  OF  WATKINS  AND  CLARY 

The  phenomenon  wherein  a  different  number  of  cir¬ 
cumferential  waves  are  present  at  the  two  cone  edges  xr>  an 
apparent  natural  mode  was  observed  by  Watkins  and  Clary 
(12)  but  not  by  Hu,  Gormley,  and  Lindholm  (data  used  in 
Chapter  V) .  As  was  mentioned  in  Chapter  I,  the  validity 
of  the  experimental  results  of  Watkins  and  Clary  was 
questioned  by  Hu  (5).  Hu  points  out  the  care  that  must  be 
exercised  in  setting  up  such  an  experiment  in  order  to 
excite  only  the  true  natural  modes  of  the  cone.  A  reason¬ 
able  theory  which  explains  the  results  observed  by  Watkins 
a  .  Clary,  based  upon  theoretical  frequencies  and  mode 
shapes,  is  presented  in  this  chapter. 

Geometry  of  the  conical  shells  taken  from  Ref.  12 
and  used  here  as  comparative  cases  are  given  by  Table  VII. 

Comparative  results  for  Cases  5,  6,  7.»  and  8  are 

shown  by  Figures  6.1  through  6.4.  The  frequency  parameter 
1/2 

A  plotted  in  these  figures  is  defined  as 
A1>/2  =  u)XL(p(l-v2)/E)1'/2 
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TABLE  VII 


CONE  GEOMETRY,  CASES  OF  WATKINS  AND  CLARY 


Case  5 

Case  6 

Case  7 

Case  8 

Cone  Half- 
Angle^ 

5SSS  SSSSSfiSSSSaStaMBSttSSSS 

3.18° 

7-59° 

14.04° 

23.96 

Major  Edge 
Radius, 

14. 0“ 

14. 0” 

14.0" 

14. 0" 

Minor  Edge 
Radius,  rQ 

12.0" 

10.0" 

8.0" 

6.0" 

Length,  L 

36.0" 

30.0” 

24.0" 

18.0" 

Thickness,  h 

0.007” 

0.007” 

0.007" 

0.007" 

Completeness 

Parameter 

s/si 

1.16 

1.40 

2.33 

1.75 

Frequency  Parameter  A 


Me.  Circumferential  Waves  m 
Figure  6.1  Frequency  Parameter,  Case  5 


m 


No.  Circumferential  Waves  m 
Figure  6.2  Frequency  Parameter,  Case  6 
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where  X  is  I/rQ.  This  is  the  parameter  used  by  Watkins  and 
Clary.  It  can  be  seen  that  the  first  mode  results  compare 
favorably  with  the  experimental  data  when  the  circumfer¬ 
ential  wave  number  observed  by  Watkins  and  Cxary  at  the 
major  (large)  end  is  used.  The  semi-empirical  equation 
proposed  by  Hu,  et_  al . ,  (6)  also  compares  favorably  to  the 
Watkins  and  Clary  data  points  associated  with  the  circum¬ 
ferential  wave  number  at  the  major  edge.  On  the  other 
hand,  the  data  which  corresponds  to  the  circumferential 
wave  number  observed  at  the  minor  (small)  end  of  the  cone 
correlate  reasonably  well  with  the  second  analytically  com¬ 
puted  mode.  This  fact  supports  the  theory  that  the  experi¬ 
mental  phenomenon  observed  by  Watkins  and  Clary,  wherein 
a  different  number  of  nodes  were  observed  at  the  two  cone 
edges,  may  have  actually  been  a  combination  of  the  natural 
first  and  second  mode  shapes  corresponding  to  different 
circumferential  wave  numbers. 

In  support  of  the  above  theory,  it  was  noted  from 
the  computed  data  that  the  predominant  displacement  of  the 
first  mode  occurs  at  the  large  end  of  the  cone,  whereas, 
the  predominant  displacement  of  the  second  mode  occurs  at 
the  opposite  edge.  Consider,  for  example,  the  cone  of 
Case  8.  Experimentally  it  was  observed  that  one  mode  had 
three  circumferential  waves  at  the  minor  edge  and  five 
waves  at  the  major  edge.  The  computed  frequency  parameter 


for  m=3,  second  mode,  is  A  7 *"  =  0.0180  and  for  m=5,  first 
1/2 

mode,  is  A  '  =  0.0195*  The  corresponding  mode  shapes  are 

presented  in  Figure  6.5*  Note  that  the  predominant  dis¬ 
placements  are  as  indicated  above.  The  closeness  of  the 
frequencies  and  the  shapes  of  the  true  natural  modes, 
therefore,  support  the  theory  that  these  two  modes  could 
easily  have  been  excited  together  in  the  experimental  case 
with  the  observed  results  that  three  circumferential  waves 
would  appear  at  the  small  end  of  the  cone  while  five  waves 
would  appear  at  the  large  end.  Many  other  similar  examples 
are  evident  from  the  data  such  as  in  Case  7  where  for  m=5, 
the  second  mode  could  have  easily  been  excited  with  the 
first  mode  for  m=7.  In  this  case,  computed  frequency 
parameters  are  0.0443  and  0.0450>  respectively.  The  cor¬ 
responding  mode  shapes  are  similar  to  those  presented  in 
Figure  6.5. 

It  can  be  shown  that  for  certain  combinations  of 
circumferential  wave  numbers,  coupling  of  the  theoretical 
modes  is  possible  as  a  result  of  the  six  point  string 
suspension  system  utilized  in  the  experiments.  This 
analysis  is  given  in  Appendix  C.  However,  only  two  spe¬ 
cific  modes  of  the  many  observed  by  Watkins  and  Clary  fall 
into  this  category.  Therefore,  coupling  due  to  the  sus¬ 
pension  is  of  minor  importance  and  does  not  alter  the 
validity  of  the  previous  discussion.  The  presence  of  the 


Normalized  Radial  Displacemen. 
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Figure  6.5  Normalized  Radial  Displacements,  Case  8 
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suspension  also  has  a  minor  numerical  effect  on  the  values 
of  the  frequencies  except  in  the  few  cases  where  the  sus¬ 
pension  and  nodal  points  coincide. 

Due  to  the  fact  that  the  experimental  frequency 
iata  plotted  in  Figures  6.1  through  6.4  correspond  to 
complex  mode  shapes,  whereas,  the  theoretical  frequencies 
represent  pure  modes  having  the  same  number  of  nodes  at 
both  edges,  comparisons  of  the  data,  other  than  those 
already  made,  are  difficult.  Since  the  mode  shapes  are  so 
drastically  different,  conclusions  regarding  a  comparison 
of  frequencies  would  net  be  entirely  valid. 
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CHAPTER  VII 


CONCLUSIONS 

The  results  of  this  investigation  have  added  cred¬ 
ibility  to  the  experimental  data  of  Hu,  Gormley,  and 
Linaholm,  and  offered  a  plausible  explanation  to  the 
phenomenon  observed  by  Watkins  and  Clary.  This  phenome¬ 
non,  wherein  different  numbers  of  nodes  were  observed  at 
the  two  cone  edges,  was  most  likely  caused  by  the  combina¬ 
tion  of  a  fundamental  mode,  having  a  given  number  of  cir¬ 
cumferential  waves,  with  a  higher  (second)  mode  having  a 
smaller  number  of  nodes.  This  is  evidenced  by  the  close¬ 
ness  of  the  respective  frequencies  and  the  location  of 
the  relative  predominant  displacements.  Other  possible 
eauses  for  the  observed  complex  mode  shapes  Include 
material  nonlinearities  such  as  might  be  caused  by  non¬ 
constant  shell  thickness  or  seams,  material  or  acoustic 
damping,  and.  In  a  few  instances,  coupling  due  to  the 
manner  of  suspension.  These,  however,  are  believed  to  be 
of  minor  importance. 

Results  of  the  higher  modes  (not  all  presented 
here)  could  be  improved  considerably  by  adding  more  terms 
to  the  assumed  polynomial  displacement  functions.  This 
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may  also  require  the  use  of  double  precision  techniques  in 
the  computer  inversion  of  the  stiffness  matrix  since  the 
accuracy  of  the  matrix  inversion  decreases  as  the  order  of 
the  matrix  is  increased.  More  sophisticated  functions 
could  also  be  used.  However,  the  integration  required  to 
determine  the  matrix  elements  may  then  become  unwieldy. 

Use  of  the  more  general  Flugge  equations  given  in 
Appendix  A  would  not  have  significantly  affected  the 
results  of  this  study  since  the  terms  omitted  in  the 
Donnell  type  equations  utilized  here  are  of  higher  order 
for  the  very  thin  shells  which  were  considered.  For  shells 
having  the  quantity  h/(s  sin  a)  greater  than  1/30,  it  is 
recommended  that  the  more  general  equations  be  used.  The 
thin  shell  assumption  h/(s  sin  a)  «  1  was  certainly  satis¬ 
fied  since  this  quantity  varied  between  1/272  and  1/2000 
for  the  eight  cases  discussed. 

Based  upon  comparisons  with  available  experimental 
data,  it  has  been  demonstrated  that  the  modified  Galerkin 
procedure  can  be  utilized  to  determine  the  frequencies 
and  mode  shapes  of  a  thin  conical  shell  with  free  edges. 

The  assumed  displacement  functions  were  found  to  ade¬ 
quately  represent  the  motion  of  the  shells  considered,  at 
least  for  the  determination  of  the  lower  modes.  A  logical 
extension  of  this  work  would  be  to  consider  the  vibrations 


of  the  conical  shell  with  other  boundary  conditions,  some 
of  which  have  not  been  investigated.  The  same  general  pro¬ 
cedure  could  be  used. 


* 


APPENDIX  A 


COMPLETE  FLUGGE  EQUATIONS 


The  equations  presented  in  this  appendix  were 
derived  from  Flugge's  (3)  fundamental  equations  for  a  con¬ 
ical  shell,  which  are  valid  for  a  variable  wall  thickness. 
His  equations  were  modified  by  1)  restricting  the  equa¬ 
tions  to  shells  having  a  constant  thickness,  2)  adding 
inertia  terms,  and  3)  changing  the  nomenclature  to  agree 
With  that  used  in  this  study.  The  equations  are  based 
upon  the  assumptions  the.t  the  displacements  are  small  and 
that  normals  to  the  middle  surface  are  preserved  as  such 
during  deformation.  The  additional  terms  in  the  following 
equations  which  are  not  found  in  Eqs.  (2.33)  to  (2.35) 
basically  represent  the  influence  of  change  in  curvature 
on  the  normal  and  shearing  forces  and  the  influence  of 
middle  surface  strains  on  the  moments.  For  very  thin 
shells  these  terms  are  of  minor  importance  and  consequently 
the  simplified  equations  were  used  in  this  study. 
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D[ 


v  esc  a  .  vu1  ,  u 


s 


w 


p 

cot  a,  _4_  _  .  Dh  r-f3-v)v//  *csc  a 
— 5 — J  cot  a  +  22“  l  J - 4 - 5 - 


2s 


(3-v)vt*csc  a  _  (3+v)v‘cbc  a  _  u^_  +  (l-v)ut ‘ *csc2a 


2s 


3 


2s 


2s- 


+  (l:yh:i^2a  .  u^cotjt,  +  u/frttqtfal]  cot  a 
2s4  s3  s4 


Dh2  r.^  .  2w"  *  *esc2a  .  w - csc^a  ,  2ww  2w'**csc2a 

LW  +  - 5 -  +  - 77 -  +  — - 5 - 


T2T 


O  O  O 

w"  .  2w*’cot  a  esc  a  ,  4w‘ ‘esc  a  .  w* 

r?  + - 3 - +  — ?r —  +3 


+  w(g-tcotatt)cotga]  .  ph  |fw  =  0 
s  dt2 


(A.3) 


kl 


FLOW  DIAGRAM:  CONICAL  SHELL  VIBRATION  COMPUTER  PROGRAM 


L=0 


4 — i  input  mm 

I  PRINT  INPUTS] 

^  _ 

I  COMPUTE  CONSTANTS! 

♦ 

i COMPUTE  SUB-MATRICES 


[ASSEMBLE  STIFFNESS1 
MATRIX  f  K] _ 

* 

I  INVERT  STIFFNESS 
MATRIX  [K] 
(SUBROUTINE) 

*  - 
COMPUTE  MASS 
MATRIX  [M] 

COMPUTE  [KJ-'lM] 
(SUBROUTINE) 

ITERATE  TO  OBTAIN" 
FREQUENCY  AND 
MODE  SHAPE  {♦} 
(SUBROUTINE) 

[COMPUTE  FREQUENCY 
PARAMETERS;  PRINT| 
OUTPUTS _ 

COMPUTE  DISPLACEMENTS] 
AND  PRINT 

tl^5  _ 

{compute  and  print 
[[♦]t[MJ[#3  for 
ACCURACY  TEST 


L  =  L  +  1 


[SWEEP  PREVIOUS 
MODES  AND  COM- 
PUTE  NEW  [U] 


COMPUTER  PROGRAM  NOMENCLATURE 


ALPHA 

ss 

a  (half  cone  angle) 

Ctf,  M 

= 

m  (no.  circumferential  waves 

DISP 

= 

displacements 

E 

= 

E  (Young’s  Modulus) 

FPAR 

1/2 

frequency  parameter  A 

FPAR2 

= 

frequency  parameter  0 

H 

= 

h  (thickness) 

L 

= 

count  on  modes 

OMEGA 

= 

frequency 

PHI 

= 

mode  shape 

RHO 

= 

P  (density) 

RL 

= 

RO 

s 

ro  1 

S 

= 

s  I 

Defined  by  Figure  2.1 

SI 

= 

S1  | 

S2 

= 

S2  i 

SMI  to  SM2  =  submatrices  used  to  generate  XK  and  XM 

VI  to  V-32  =  various  combinations  of  s,,  sp,  v,  a,  etc. 
which  are  used  repeatedly1  d 

XK  -  stiffness  of  matrix 

XKI  =  inverse  of  stiffness  matrix 

XLEN  =  L( length) 

XM  =  mass  matrix 

XMU  =  v  (Poisson’s  ratio) 


,  <Kt»,  /<  /  ‘  »*4H>  V  *>'«’<■*  >  r  if**  I  ,  »  «  *«>*,  „  .  „  »  •.  ,.  ^  ,,\  /  r  ,t,«  1^,  >t  ,  ,<  » J,  *"  K  ‘-,1  *i  •  fr  -  **$'*&•> 
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C  THE  PURPOSE  OF  .  THIiLJ’RQGRAil--iJ>  IM£_-SQJJtfICN_OF_HULr.I-OEGR£t 
C  OF  FREEDOM  VIBRATIONS  PROBLEMS  FOR  THE  NATURAL  FREQUENCIES 

C  AND  MODE  SHAPES  BY  ITERATIVE  METHODS. 

C  BASIC  EQUATIONS  U/(GMEGA**?1*A«XK1*M*A 

C  WHERE  A  ON  THE  RIGHT  IS  A  TRIAL  EIGENVECTOR 

C  XKI  IS  THE  INVERSE  OF  THE  STIFFNESS  MATRIX 

C  M  IS  THE  MASS  MAIfUX  _  _  _ _ 

C  A  ON  THE  LEFT  IS  THE  UNITIZED  RESULT 

C  SIZES  XK=M*M  XM=M*M  A*M*I 

0 1 ME  N  S I ON  XK 121*21) .XKI (21 .21 1 .UI21.21 1 «UO!21,2l I ,XM( 21,211 ,C 1 21* I ) 
It  PHI 1 21*21) .PHI T{  21*21) » SMI I A #4) ,SM2(4.4) ,SM3t4,4| *SM4(4*4)« 

?SM5( 4*4) «  SM61 4*4) • SM7 14*4) • SM8 (4 .4) *SM9 (4 *4 ) *SM10 14* 4) • SMI K 4,4)  \ 

DIMENSION  RO«121*211cQLSPt2JL*21) .  _  . 

100  FORMAT! 1H1*4X«18HENTER  NEW  DATA  SET,//) 

101  FORMAT! 12) 

102  FORMAT! 10X.2HM-, 12) 

103  FORMAT! 8F 10. 0  ) 

10?  FORMAT! lOXfllHMOOE  SHAPE  ,I2.5X,2H0  ,I2,IPE15.7) 

108  FORMAT!  10X,6HQMEGA.a.I2*ll?E2Q. 7)  .  . . _  . 

109  FORMAT! 5X.14HEND  OF  PROBLEM) 

110  FORMAT!//) 

111  FORMAT  1 10X.26H0RTHQGCNALIZED  MASS  MATRIX,2I2,1PE15.7) 

113  FORMAT! 10X. 19HFRE0UENCY  PARAMETER. 1PE15. 7) 

114  FORMAT! 1H0.10X.44HFREQUENCIES  AND  MODE  SHAPES  OF  CONICAL  SHELL) 

115  FORMAT!  IX, IP 4E1Z. XI _ _  _  _  _ ... 

116  FORMAT! 10X,6HLENGTH,l0X,9HTHtCKNESS.8X,9HNO.  WAVES. 12X. 1HE ) 

117  FORMAT! 10X. 2HS1 • 15X,2HS2. 13X*5HALPHA) 

118  FORMAT! 10X.2HR0, 15X,2HRL.15X.3HRH0.14X,3HXMU) 

120  FORMAT! 1QX.4HS/S2, 14X ,1HU, 15X ,1HV, 15X ,IHW ) 

C  READ  INPUT  DATA  AND  COMPUTE  XKI 

C  INPUT  DATA _  ... _  _ _  _ _  ...  .  _... 

PRINT  114 
l  READ  101, M 
PRINT  100 
L»0 
NOTE® 1 

PRINT  1Q2.M. _ _ _ _ _ _  _ _  ...  . 

PRINT  110 

2  READ  103,R0,RL.RHO,XMU,XLEN,H*CW,E 
PRINT  110 
PI=3. 1415927 

ALPHA=ATANFI IRL-RGl/XLENI 

S1*R.0/SINE1ALEHA1 _ _  _ _ 

S2«RL/SINF! ALPHA) 

ALP0*ALPHA*57. 29578 
XLAM®XLEN/RO 
PRINT  113 

PRINT  115.R0,RL,RH0,XMU 

PRINT  life _ ... _ _ _  ...  _  .. 

PRINT  11 5,XLEN,H,CW*E 
PRINT  117 

PRINT  115,S1,S2.ALPD  j 

PRINT  110  .  .  1 

XT=!H*H1/12.0  1 

_ £  C.QMEU  INCONSTANT  VALUES _ _ _ _  1 
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VJ*S2-S1  ,  - 

V2sS2*S2-Sl*St 

V3=S2**3-S1**3 

V«=S2**4-S1**4 

V5*S2**5-S1**5 

V6*S?**6— Sl**.6 

V7*LQGFC  S21-LOGF.I  Sll  .  .  _ _ 

V8*S2*S2 
V9«S2**3 
V 10* S 2*44 
VI l*CR/SINFI ALPHA) 

V12*V1!*V11 

V13*C0SF<ALPHA)/SINFIMPHAI  .  . 
V14*V13*V13 

Vl5*l.0*0.5*( 1.0-XMU) *V12 
V16*Q.5*< 1 • C—XMU) ♦ V12 
V17*(XT*V12I*I V12-4.0) 

Vi 8*— ( 1.0-XMU! /2.0 
V19*XT*(2.C*Vl2U..fl)  _ 
V2Q=-XT*V12*(3.0-XMU» 

V21*-XT*XMU*V12 
V22*XT*XMU 
V23=1.0/S2-l.0/Sl 
V 24*1 .0/1 S2*S2 )-l. 0/ ( Si*Sl ) 
V25*VU*(3.0-XMU)/2.XL  .  . 

V26=— Vll*( 1.0+XMU) /4.0 

V27=-Vll*U.C-3.  0*XMUI/4.0 

V28*  (  RHO*  <  1.0-XMU*XMU) 1/16*386.088) 

V29*S2**5 

V30*S2**6 

V3l*S2**7-Sl**7 _  ...  _ ... 

V32=S2**6-S1**8 
C  COMPUTE  SUBMA TRICES 
N=M/3 

SM1< 1.1) *V2/2. 0 
SM1C 1,2)*V3/C3.0*S2) 

SHM  1*3)*V4/(4.0»V8) _ _ 

SM1( 2.21= SM 1(1*31 
SMl(2,3J-V5/(5.n*V9)  _ 
SMI(3*3I=V6/(6.0*V10I 
SMI I 1,4)*SM1<2*3) 

SM 1(2*41* SM 1(3*3) 

SHU  3* 41=  ¥11  _ 

SMI I  4*4 ) =V32/( 8. 0*  V30) 

DO  20C  l*l*N 
DO  200 

200  SM1I J« I )*SM1(I *J) 

SM2C 1.1)*V7 

SM2< lt2)*Vi/S2  _ _ _ _ 

SM2I i.3)=V2/I2.0*V8) 

SM2I 2«2)*SM2( l«  3) 

SM2(2.3)*V3/(3.0*VS) 

SM2I 3,3)*V4/(4.0*V10) 

SM2I 1*4I*SK2(2*31 

SM2(2.4)*SM2<3.3) _ 

SM2I 3»4I*V5/(5.0*V29) 


'•  K» 
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OQ  210  1*1. N 
00  210  J*I *N 
210  SM2‘  J.H=SM2(I.J) 

SM3( 1.15=0.0 

sr.ii  i .2)  =c. o 

.SMJ1X^U*£U0 _ 

SM1I 2,2)*V2/(2.0*V8) 
SM3(2,2)=(2.  Q*V3l/(3.  Q*.V9I 
SM3(3*3)=V4/V10 
SM3(l.Ar-0*0 

SM3<2.4I=<3,0*V4)/(4.0*V10» 
SM3L3j4.1_*ii>.Q*V5)  /  (5 .0*1/29  L 
SM3U.4)  *(3.0*V6)/(2.0*V30) 
00  220  1*1, N 
PO  220  J=I.N 
220  SM3! J. I 1=SM3I I ,J1 
SM4< l.ll=-V24/2.0 

SM4U,21~^V23/-S2 _ : _ 

SM4< 1*3)=V7/V8 
SM4( 2,2 )*SM4ti«3J 
SM4{ 2,31 *V1/V9 
SM4I 3,3t*V2/(2.0*V101 
SM4I 1,4)*SM4<2,3) 

SM4( 2,41 =SM4(3«31 _ 

SM4I 3.4) *V3/( 3.0*V29) 
SM4C4,4)*V4/I4.G*V30) 

00  230  1*1, N 
DO  230  .1*1, N 
230  SM4I J,I»*SM4(I,J) 

00-240-  L»UN _ 

00  240  J*1.N 
240  SM5( I«J)*0.0 

SM5I3, 31*12. C*V2)/V10 
SM5{  3,41  *14. G*V3) /V29  ..  .. 

SM51 4»3)*SM5f3,4) 

.  _  SM-§1Al,41.*( 9. Q* V41/V3Q _ 

SM6I 1,1)*0.0 

SM6t 1«2)*V1/S2  . 

SM6( 1»3I*V2/V8 

SM6(  2,21*0.0  -  ...... 

SM6( 2,3)*V3/(3.0*V9I 

_ .  Jl _ 

SM6C 1,  V3/V9 

SM6(2f4)*V4/J2.0*yi01  .  „  . 

SM6I 3,41 *V5/( 5.0*V29f 
SM6I4, 41*0.0 
00  25C  1*1, N 


250  SM6( J , I )*— SM6( I , J) 
SM7( 1,11*0.0 
SM7C 1.21*0.0 
SM7< 1,31*0.0 
SM7(2»2)*V7/V8 
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SM74 1,  4)=0._Q _ 

SM742.4)=<3.C*V2)/42.0*V10) 
SM74 3.41=1 2. C*V3)/V29 
SM744. 41*49. 0*V41/44. C*V30) 
DO  260  I =1 «N 
00  260  J  =  I.N 

260  SMZ4  j  .11  =SKXLi,JJ _ 

SM84 1.11=0.0 
SM8( 1«21*V1/S2 
SM8( 1. 3I*V2/V8 
SM8( 2.21= SM 8(1.3) 

SM84 2*  3 1=V3/V9 

SMS  (3,  31*Y4mQ _ 

SM84 1.4  1*SM842«3) 

SM8I 2.4)=SM8(3«3) 

SM84  3.41*V5/V29 
SM8(4.4)=V6/V30 
DO  270  1=1. N 

.  OC.  27C  .  O=l*N _ 

270  SM8U.I 1=SM84I .J) 

SM94 1.1 )*V24 
SM94 1.21=V?3/S2 
SM94 1.31=0.0 
SM94  2.21=0.0 

-  _ 

SM94  3. 3I=V2/V10 
SM94 1.41»SM942.3) 

SM94  2.41=SM943.3) 

SM94  3»41=V3/V29 
SM94  4.41*V4/V30 

. PQ-2aO_l=  1..N _ 

DO  280  J=I.N 
280  SM94  J,n=SM941  .J1 
SM104 1. 1 1=0.0 
SM104 1. 2)=V23/S2 
SMI04 1.31=0.0 

JjMlQ12*2JL?.Q..JS _ 

SM104 2.31=4 3. 0*Vll/V9 
S  M 1 0  4  3  •  3 1  *  4  4  .  0*  V2 1  J.  yi  0. . 

SM104 1,41*4 3.0*Vli/V9 
SM104  2t41*44.04y2i/.ViO.. 

SM 10 4 3,41*4  5.04V31 /V29 
.  SM1Q44. 41=46, Q4V41/V30 

DO  2 90  (*1.N 
00  290  J»I,N 
290  SK104J.ll»SM104I.Jl 
00.  300  1  =  1, N 
00  300  J=1.N 

300  SMll 4  I. J 1=0. 0 _ 

SM1142, 31=42. C4V1I/V9 
SMI 14  3»21*SM 1142,31 
SMI 14 3.31=4 4.0* V2I /VI 0 
SM  11(2. 41  *4  3.06V2J./V10 
SMI 14 4.21 «SM1 14 2. 41 

_ SM1143, 41  =  46. Q4V31/V29 

SM1144. 3)=SM1143.4) 
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SMULi4.4i  =  l/i.0*V41  A3C _ 

COMPUTE  STIFFNESS  MATRIX 
ALPHA**  ALPHA  SUBMATRIX 
DO  310  1=1. N 
00  31 C  J* l.N 

XKI  I,  JI«V15*SM2U.  JI»SM3(I.J)+XMU*SM8II.JI 

B£IA=&£1A  SUBMATRIX _ 

IPN'  I*N 
JPN*J+N 

XKI IPN. JPN  J =V16*SM2I I . J»-V18*I SM3 If . J1-SM8I I • J)  ) 

GAMMA-GAMMA  SUBMATRIX 

IP2N*I«-2*N 

JP2tW*2*N _ 

XKI IP2N© JP2NI*V14*SM2! I  .  J1  +  V17PSM4II .  J) 
im9*SM7ll.Jl>XT*SM5IIt.Jl*V20*SM9II.J)*V21*SHlOU., 
2+V22*SMl If  I  * J  t 
ALPHA-BETA  SUBMATRIX 

XKI I . JPN)=V25*SM2 ( I . J»*V26*SM6 1 1 • J) ♦V27*SM8 ( I.J) 


XKII. JP2NI*V13*C-SM2I I  .Jl  +  1  XMU/2.0I*<S*6U  .  J1-SM81 
BETA-GAMMA  SUBMATRXX 

XKI  f.PN.JP2N)=-VU*V13*SM2U.J> 

00  3 20  I*1.M 
on  32C  J*I «M 

XKIJ. I1  gXKJIt.il _ _ _ 

CALL  MAINV! XK.XKI «M) 

CALL  MAMUL! XKI .XK.U.M.M.M) 

COMPUTE  MASS  MATRIX 
00  350  I*l*N 
00  35C  J*1.N 

I£i4*LfN _ 

JPN*J+N 

JP2N=J*2*N 

IP2N=l*2*N 

XMLI. Jl*V28*SMit  l.JJ,  .  ... 

XMI IPN. JPN )*XMI I.J) 


i. jin 


COMPUTE  U*XK!«XM 
CALL  MAMUL€XKi*XM.U«.M»M.M> 
ITERATE 

5  CALL  ITERIU.O.OMEGA.M.MOTE) 
L*L*1 


PHI(I.L1*0(I 1 

10  PHITIL.I  )*0U  1  . . . 

FPAR*SQRTFI 1 1 1.0-XMU*XMUI •RHOPXLAMPXLAM*XL£NPXLEN*OMEGAPQMEGA) 

l/(£*386»088t| . 

FPAR2«I!FPAR*RL*R01/1XLEN*XLEN» 


PRINT  110 
PRINT  10Q.L. OMEGA 
PRINT  1 13.FPAR 
PRINT  11 3.FPAR2 
PRINT  110 


00  60  J*1.3 
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m 

Sr 


-£ 

2 


60  RQWU,I1*G. . ..  .  ..  .  . 

PRINT  120 
FL*Sl/$2 

SINC*< S2-S1 1 /( 20*06 S2 ) 

6?  ROWU, 11*1.0 
ROW! i,21*FL 

ROWI  1*31 =FL*Ei _ _ _ ._  .  . _ _  _ 

ROWI  1.41«FL*FL*FL 
ROM! 2,51*1.0 
ROWI 2.61-FL 
ROW! 2,71*ROMI 1.31 
ROWI 2.81 aROW( 1.4) 

ROWI  3*91  rl«Q  .  . .  ..  ...  . 

ROWI 3*10) *FL 
ROWI 3*11 l=RQWf 1*3) 

ROWI 3«  12 1 =ROWl 1*41 

CALL  MAMULI RGW.Q.DI SP.3.M.11 

PRINT  ll5.FL.OISPIl.ll.DISPI2.il •0ISPI3.1) 

FL*FLfSlNC _  _  _ .... 

IFIFL-1.0162,63,63 
63  CONTINUE 
PRINT  110 
1FIL-51  6.8,8 
C.  FORM  NEW  U=U-U0 

6  CALL  MAMULtXH«0«.Uftj.M.*.Mtll  ....  ...  . 

CALL  MTMULIO.UQ.XK.l.M.ll 
DEN*XKIl.ll*|OMEGAI**2 

CALL  MTMUL  (0,XM,XK.1,M.M1 
CALL  MAMULIG.XK.UQ.P.l.Ml 
00  7  1*1. M 

DO  1  4*1,  M _  _  .  . . 

UOU  ,JI=UQ(  I.JI/DEN 

7  Ult.JIMJII.JI-UQU.JI 
AA*0. 

BB»C. 

00  11  1*1. M 

AA*AAU, _ _  _  _  _ _ 

11  BB*BB+Q( 1 1 

IF  IAA-BBI  5,12,5 

12  NQTE*2 
GO  TO  5 

C  FIND  ORTHOGONAL I ZEO  MASS  MATRIX  TO  TEST  FOR  ACCURACY 

8  CALL  MAMUL_I  XM.PHI  ,XK,M_>H.5) _ _  _  _ 

CALL  MAMUL  IPH1T.XK.XKI .5.M.51 

PRINT  lll.(II.J,XKI(I.J).I*l«51«J*1.5J 
PRINT  110 
PRINT  109 
GO  TO  l 
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SUBROUTINE  MAM'JL  IA*B*C*M*N,MM» 
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SUBROUTINE  i!Ut!!UL-JAtilfc£ti4*NxRJU _  _ 

THIS  SUBROUTINE  IS  FOR  THE  COMPUTATION  OF  MATRIX  C  FROM  C«A*B 
WHERE  CIS  AN  M*MM  MATRIX*  AIS  AN  M*N*  BIS  AN  N*MM 
DIMENSION  A I ?l *21) *8(21*21) tC<21»21) 

DO  3  1=1. M 
DO  3  J=ltMM 

CII.Ji*0-  _ _ _ _ _ _  _ _ 

DO  3  K*l«N 

3  C(I*J)=C(I»J)+AII.K)*B(K.J) 

RETURN 

END! 1.0«0*0«0*C*l*Ct0*l*0*0*0*0*0) 


wr>)<  AwHW*Ml  ►*»«»•»■  <i<l»  W»s«- 


SUBROUTINE  MAINV  (A.AINV.M1 
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%  SUBROUTINE  MALNV._IJUA1NV.ML  _  . 

C  THIS  SUBROUTINE  IS  FOR  COMPUTING  MATRIX  AINV  WHICH  IS  THE 
C  INVERSE  OF  MATRIX  A 

••  '  DIMENSION  A(21,21l  .8(21,421  .AINVI21.21) 

N  *  2*M 

*  DO  4  I  *  1 ,  M 

DO  5  J=1.M  _ _  _ _ _ _  ..  .  - 

B<  I.  J)  *  A  ( I  «  J ) 

5  CONTINUE 
4  CONTINUE 

L*M+1 

X  DO  6  1*1, M 

DO  7  J»UN  . . . .  ...  .  .. 

IF  II-J+MI  8,9,8 

8  B (  I, J|  *  0, 

J  GO  TO  7 

9  B(  I,  J  )  *  i, 

7  CONTINUE 

6  CONTINUE . .  .  _  ...  . . . . 

*  DO  10  J=l.M 

C  *  ft( J, J  ) 

IF  (Cl  20.21,20 
20  00  11  K*  l  «N 
11  BIJ.K)  *  BIJ.KI  /  C 

5  DO. 14  L=l.M _  _  . . . 

3  IFIL-JI  13,14.13 

s  13  D  *  BIL.JI 

{  *  DO  15  K*l,N 

|  15  BIL.K)  *  B(L*K)  -  BIJ.KI  *  D 

I  14  CONTINUE 

!  10  CONTINUE. _ _  _  _ _ _ 

*  L  *  M*1 

i  DO  16  1*1. M 

l  DO  17  J*l,M 

L*M+J 

f  A  INV( I.J)*B(I,L) 

17  CONTINUE _ _ 

16  CONTINUE 
I  GO  TO  23 

■  21  PRINT  22 

22  F0RMATI///.5X.39HINVAUD  SOLUTION — ZERO  DIVISOR  IN  MAINV,///) 
)  23  RETURN 


L 
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SUBROUTINE  ITER  1U.Q. OMEGA, M. NOTE!  82 

SUBROUTINE  ITER  Hi. Q* OMEGA, M.NQIE) 

DIMENSION  U< 21,2 II  *0(21,1)  ,A(2l,l> 

X*0. 

on  i  i*i, m 

l  om*i. 

GO  TO  I  3, 7), NOTE 
7  DO  8  1=1«M«2 

a  cm=-i. 

3  CALL  MAMUL(U,0,A,M,M,1) 

RQ*A(M! 

8*0. 

00  6  1=1. M 

6  8=fiKUI) 

C=0. 

DO  4  1=1, M 
0< I }=A( I i/RO 

4  c=c*om 

X=X*l. 

IFIX-100.1  10*10.5 
10  IF  lABSFlB-Cl-. 0000001)  5.5,3 
5  PRINT  300. X 

IF{ R0  111.  12.12 

11  PRINT  301.RO 

12  RO=ABSF!ROI 
QMEGA=S0RTF(  1«JQ/R0I 

300  FORMAT  ( 10X. 21HNUMBER  OF  ITERATIONS  .1PE15.7.//1 

301  FORMAT! 10X, 14H0MEGA  NEGATI VE.1PE15.7,/) 

RETURN 

END! 1.0.0,0.0,C,1,C. 0,1 ,0,0, 0,0,0! 
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APPENDIX  C 

COUPLING  OP  MODES 

A  mathematical  demonstration  that  the  natural  modes 
of  the  conical  shell  are  coupled  by  the  six-point  string 
suspension  utilized  by1  Watkins  and  Clary  is  presented  in 
this  appendix. 

In  the  expansion  of  the  kinetic  energy  expression 

T  =  1/2  JJph(u2+v2+w2)s  sin  a  dsd8  (C.l) 

the  cross  terms  involving  m  and  n,  two  distinct  modes,  inte¬ 
grate  to  zero,  verifying  that  no  inertial  coupling  exists. 
Elastic  coupling  due  to  the  presence  of  the  strings  will  now 
be  examined.  Consider  the  six  support  strings  to  be  located 
at  8=0,  tt/3,  2tt/3,  4tt/3  and  5tt/3  •  The  displacement  of  the 

4*  Vv 

strings  in  the  nr  mode  is  longitudinal  only  and  is  given 

by 

U1  "  %  qi  =  am(s)  cos  ra9’ql  (c-2) 

m  " 

u2  =  q2  =  am(s)  sin  (c-3) 

The  two  functions  are  needed  for  complete  generality  in 
orientation  of  the  mode  with  respect  to  the  support  strings. 
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The  total  displacement  In  the  direction  of  the  suspension 


strings  is  then 


u  =  +  Ug 


(0.4) 


Forces  in  the  strings  are  given  by  Eqs.  (C.5)  to  (C.10) 
where  k  represents  here  the  string  spring  constant. 

6  =  0:  -  k^  -  ku2  =  -  kam(s)[cos  0*q1  +  sin  0-q2]  (C-5) 

0  =  ir/3 :  -  ku-,  -  ku2  =  -  kam(s)[cos  m7r/3 •  q]_  +  sin  imr/3*q2] 

(0.6) 

9  =  2tt/3:  -  kux  -  kti2  =  -  kam(s)[cos  2mir/3-q1  +  sin  2imr/3'q2] 

(C.7) 

0  =  ir:  -  ku1  -  ku2  =  -  kam(s)[cos  mir^  +  sin  nnr-q^ 

(C.8) 

0  =  4ir/3:  -  ki^  -  ku2  =  -  kam(s)[cos  4m?r/3 •  +  sin  4imr/3*q2] 

(0.9) 

0  =  5tt/3:  -  ku2  -  ku£  =  -  kam(s)[cos  5mir/3 •  q±  +  sin  5mTr/3'q2l 

(C.10) 

A  virtual  displacement  in  the  nth  mode  is  given  by  Eqs. 

(C.ll)  to  (C.16). 

0-0:  6u^  +  6u2  an(s)[cos  O-Sq^  +  sin  0'6q^j  (C.ll) 

6  =  tt/3:  6^  +  6u2  =  an(s)[cos  wr/3’6q^  +  sin  rar/3*6q^] 

(C.12) 
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8  -  2tt/3:  6u^  +  6u2  =  an(s)[cos  2ror/3*6q3  +  sin  2nir/3*6q .^] 

(C.13) 

8  =  it:  5Uj  +  6u2  =  an(s)[cos  n7T*6q3  +  sin  n-(r-6q^]  (C.l4) 

8  =  4ir/ 3:  6u1  +  6u2  =  an(s) [cos  4ror/3*6q3  +  sin  4nir/3*6q^] 

(C.15) 

8  =  5ir/3:  6u1  +  6u2  -  an(s)[cos  5n7r/3-6q3  +  sin  5nir/3.6q^] 

(C.l6) 

The  virtual  work  done  by  the  string  forces,  Eqs.  (C.5)  to 
(C.10),  during  the  virtual  displacement,  Eqs.  (C.ll)  to 
(C.l6)  is  given  by 

6W  =  -  kam(s)an(s)[(qi  +  0*q2)(6q3  +  O-Sq^) 

+  (cos  imr/3 •  qx  +  sin  m7r/3-q2)  (cos  rar/3*6q3  +  sin  mr/S-Sq^) 

+  (cos  2m7r/3-q1  +  sin  2nnr/3-q2)  (cos  2mr/3-6q3  +  sin  2mr/3-6q^) 
+  (cos  m-jr-q-j^  +  sin  m7r*q2)(cos  nrr*6q3  +  sin  mr-6q^) 

+  (cos  4m7r/3  •  qx  +  sin  4m7r/3*q2)  (cos  4rar/3  ♦  6 qq  +  sin  4mr/3*6q^) 
+  (cos  5ni7r/3,q1  +  sin  5nnr/3 -q2)  (cos  5mr/3'6q3  +  sin  5mr/3.5q4) 


(C.17) 
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which  can  be  put  in  the  form 

6W  =  -  kam(s)an(s)[a1q16q3  +  a^Sq^  +  a3q2&q3  +  a^qgfiq^ 


+  higher  order  terms]  (C.l8) 

Evaluation  of  the  coefficients  a a2,  a3,  and  a^  in  Eq. 
(C.l8)  reveals  that  for  certain  combinations  of  m  and  n, 
they  are  not  zero,  indicating  that  coupling  exists.  The 
following  modes,  in  addition  to  others  involving  higher  m 
and  n,  are  coupled:  m=2,  n=4;  m=>-2,  n=8;  m=4,  n=8;  m=5,  n=7> 
m=5,  n=llj  m=7i  n=ll;  m=8,  n=10.  Since  the  functions  in 
the  coefficients  repeat,  there  are  infinite  numbers  of 
coupled  modes. 


SYMBOLS 


D 

E 

EjJ  &£>  *  *  * 

Fu,  Fv,  Pw 

h 

K 

k 

L 

M 

m 

M 


Rigidity  parameter  defined  by  Eq.  (2,l6) 

Young’s  modulus  of  elasticity 

Error  functions 

External  forces 

Thickness 

Stiffness 

Thickness  parameter  defined  by  Eq.„(2.64) 

Length 

Mass 

Number  of  circumferential  waves 
Generalized  boundary  force 


Ns,  Qs,  Mg, 


Stress  resultants 
Generalized  boundary  force 
Cone  radius  at  major  end 
Cone  radius  at  minor  end 
Generalized  boundary  force 
Longitudinal  coordinate 
Generalized  boundary  force 
Time 

Components  of  displacement 
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U,  V,  w 
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Fourier  coefficients 
Strain  energy 
Work 

Cone  semi-vertex  angle 

Function  of  s  in  assumed  u 

Function  of  s  in  assumed  v 

Frequency  parameter  defined  by  Eq.  (5.1) 

Function  of  s  in  assumed  w 

Shear  strain 

Error  functions  of  s  and  e 

Normal  strain 

Circumferential  coordinate 

Curvature 

Poisson’s  ratio 

Density 

Frequency 

Frequency  parameter  defined  by  Eq.  (4.1) 
First  partial  derivative  with  respect  to  s 
First  partial  derivative  with  respect  to  0 
First  derivative  with  respect  to  time 
Variation  of  (  ) 

Higher  partial  derivatives  with  respect  to  s 
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be  uncUasified.  Each  paragraph  of  the  abstract  shai'  end  with 
an  indication  of  the  military  sapurity  classification  ot  the  in¬ 
form?  yon  in  the  paragraph,  represented  as  (TS).  (S).  (C).  or  (V) 

There  it  no  limitation  on  the  length  of  the  abstract.  How¬ 
ever,  the  suggested  length  is  from  150  to  225  words. 

14.  KEY  VORDS:  Key  worde  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.  Key  words  must  be 
selected  to  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con¬ 
text.  The  assignment  of  links,  rules,  and  weights  is  optional.  | 
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